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£SJ ■ Abstract 

The recent introduction of a boundary stress tensor for asymptotically flat spacetimes enabled a new con- 
struction of energy, momentum, and Lorentz charges. These charges are known to generate the asymptotic 
symmetries of the theory, but their explicit formulas are not identical to previous constructions in the liter- 
ature. This paper corrects an earlier comparison with other approaches, including terms in the definition of 
the stress tensor charges that were previously overlooked. We show that these terms either vanish identically 
S^ ' (for d > 4) or take a form that does not contribute to the conserved charges (for d = 4). This verifies the 

earlier claim that boundary stress tensor methods for asymptotically flat spacetimes yield the same conserved 
charges as other approaches. We also derive some additional connections between the boundary stress tensor 
and the electric part of the Weyl tensor. 
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1 Introduction 

Conserved quantities play an important role in our description of physical systems. Particularly interesting 
cases arise in gauge theories, where conserved charges may be expressed as surface integrals at infinity. For a 
theory that is diffeomorphism invariant, like general relativity, the construction of conserved quantities depends 
crucially on the asymptotic structure of spacetime. In the case of asymptotically flat spacetimes in d > 4 
dimensions, one is most interested in conserved charges associated with the generators of asymptotic Poincare 
transformations. 

There are many different methods for constructing conserved charges in asymptotically flat gravity [1, 2, 3, 
4, 5, 6, 7, 8, 9, 10]. The various approaches offer different perspectives, but ultimately lead to equivalent results. 
For example, Arnowitt, Deser, and Misner (ADM) [1, 2] based their construction on the initial value formulation 
of general relativity. This introduced a notion of asymptotic flatness at spatial infinity by demanding that the 
Cauchy data on some initial surface (rather than the geometry of the full spacetime) should asymptotically 
approach that of a corresponding surface in Minkowski space. In contrast, a fully covariant approach to 



asymptotically flat conserved quantities was introduced by Ashtekar and Hansen in the late 1970s. Several 
important results were established within this framework, including the introduction of a precise notion of 
asymptotic flatness at spatial infinity [3, 4], and showing that the ADM 4-momentum is obtained from the 
past limit of Bondi 4-momentum [11]. More recently, a definition of conserved charges was introduced that 
makes use of a boundary stress tensor [10]. Boundary stress tensor methods for asymptotically flat space were 
motivated by the success of such methods in asymptotically anti-de Sitter space (see [12, 13, 14] as well as [9] 
and many succeeding works), which were in turn inspired by the methods of [15]. Earlier steps toward using 
such methods in asymptotically flat space were taken in [9, 16, 17, 18]. 

The starting point of the boundary stress tensor method is a valid variational principle for asymptotically 
flat gravity. The standard form of the action, comprising the Einstein-Hilbert and Gibbons-Hawking- York 
terms, must be modified because it is not stationary with respect to arbitrary field variations that preserve the 
boundary conditions. This is directly related to the problem investigated by Regge and Teitelboim in [2], who 
found that the Hamiltonian formulation of general relativity with asymptotically flat boundary conditions is 
only well-defined if it includes the additional surface term first proposed in [1]. The same conclusion applies 
to the Lagrangian formulation of the theory, but the necessary boundary terms were not investigated until 
recently; see [10] for a discussion in the 2nd order formalism or [19, 20, 21] for a discussion in the first order 
formalism. Given a valid variational principle, the boundary stress tensor T a ^ is defined by varying the on-shell 
action with respect to the boundary metric, as in [15, 12, 13]. Any asymptotic Killing field £ a can then be used 
to define a current T a & £ , and the flux of this current across a cut C of spatial infinity is a conserved charge. 
According to the results of [22] these charges necessarily generate the asymptotic symmetries of the theory. In 
[10] it was argued - on general grounds - that the conserved quantities defined in this way must agree with 
other definitions. To illustrate this point the boundary stress tensor definition of energy and momentum at 
spatial infinity were shown to agree with the charges defined by Ashtekar and Hansen in [3]. 

The comparison between the boundary stress tensor approach of [10] and other methods in the literature was 
extended in [23], where two important results were established. First, it was shown that the canonical (space + 
time) reduction of the covariant action yields the ADM result [1] with the appropriate boundary terms. Second, 
using the form 

Tab = -^—77 (nab - ftab) (1-1) 
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of the boundary stress tensor presented in [10] , explicit formulas for boundary stress tensor charges associated 
with asymptotic translations and Lorentz transformations were shown to agree with those given by Ashtekar 
and Hansen in [3]. In (1.1) Tr a b is the conjugate momentum of the gravitational field and 7r a 6 is an analogous 
contribution from the counterterms. However, it turns out that there are several additional terms (beyond those 
in (1.1)) which potentially contribute to T a b but were overlooked in [10] 1 . 

The primary purpose of this paper is to address the terms in the stress tensor that were overlooked in 
[10, 23]. Specifically, we want to resolve the tension between these 'extra terms' and the observation in [23] 
that the stress tensor (1.1) leads to charges that agree with the definition of Ashtekar and Hansen [3]. To this 
end, we explicitly calculate the extra terms and show that they do not contribute to the conserved charges. 
In d > 4 spacetime dimensions this is because the extra terms are identically zero. However, for d = 4 the 
situation is more subtle; the extra terms are not zero, but they do not contribute to the conserved charges. We 
also establish several additional results that we summarize here in simple physical terms. 

• For d > 4 the boundary stress tensor agrees with the electric part of the Weyl tensor at both first and 
second order in an asymptotic expansion. 

• For d = 4 the first order term in the asymptotic expansions of the boundary stress tensor and the electric 
part of the Weyl tensor are the same. This result also applies to spacetimes carrying non-zero NUT 
charge. At second order the boundary stress tensor is given by the second order term in the electric part 
of the Weyl tensor, plus contributions from the extra terms described above. The extra terms vanish for 



1 We thank Julian Le Witt and Simon Ross for bringing this fact to our attention. 



some spacetimes and are non-zero for others. However, in no case do they contribute to the conserved 
quantities. 

• For d = 4, with zero NUT charge, there is a non-trivial identity relating terms in the asymptotic expansions 
of the electric and magnetic parts of the Weyl tensor 

*ffi " ° E % = Cc d (aD C P4 . (1.2) 

This identity helps one see more clearly the equivalence between the boundary stress tensor and Ashtekar- 
Hansen charges. It can also be used to simplify the analysis of [23]. 

The rest of the paper is organized as follows. We begin with various definitions and a brief review of 
asymptotic flatness in section 2. We then present the covariant counterterm action of [10] in section 3, and 
vary it with respect to the boundary metric to obtain the full boundary stress tensor. In sections 4.1 and 4.2 
we discuss the asymptotic expansion of the stress tensor in d = 4 and d > 4 dimensions, respectively. These 
sections refer to calculations that are presented in detail in B and C. The stress-tensors are then used to 
construct conserved quantities at spatial infinity that agree with the usual definitions [1, 2, 3, 4, 5] of energy, 
momentum, angular momentum, and boosts. We establish this by explicitly showing that the extra terms 
overlooked in [10, 23] do not contribute to the conserved charges. Finally we consider some examples in section 
5 and close with a brief discussion in section 6. Throughout we use the notation of [10, 23, 24]. 

2 Preliminaries 

This section sets the stage for our later work by providing relevant definitions and review. In the first part 
we state our definition of asymptotic flatness in four and higher dimensions, followed by a discussion of the 
asymptotic field equations. The calculations in this part follow the approach of Beig and Schmidt [25, 26]. The 
second part contains definitions and a few useful results concerning the electric part of the Weyl tensor. 

2.1 Asymptotic Flatness 

We begin by reviewing the coordinate-based definition of asymptotic flatness that was used in [10]. This 
approach is inspired by the covariant phase space treatement of [27] and is modeled on the definitions given 
in [25, 26]. Our definition is particularly close to that of [5], which treats spatial infinity as the unit time-like 
hyperboloid. 

A spacetime of dimension d > 4 is asymptotically flat at spatial infinity if the line element admits an 
expansion of the form 

2cr 



ds 2 = ( 1 + 



p d-3 



+ 0(p-( d -V)) dp 2 + p 2 (h<2 + ^3 + 0(p-( d -V)\ dr, a dr, b +p {0{p-^-^)) a dpdrf, (2.1) 



for large positive p, where h°} and rf are a metric and the associated coordinates on the unit (d — 2, 1) 
hyperboloid H, and a, h^ are respectively a smooth function and a smooth tensor field on 7i. The coordinate 
p is the hyperbolic "radial" function associated with some asymptotically Minkowski coordinates x a through 
p 2 = Vab xa x b - In (2.1), the symbols 0(p~i d ~ 2 )) refer to terms that fall-off at least as fast as p~i d ~ 2 ) as one 
approaches spacelike infinity, i.e., p — > +oo with the coordinates rf fixed. 

Following [25, 26] we partially fix the gauge and bring the metric into the form 

ds 2 = N 2 dp 2 + h ab dri a dri b (2.2) 
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dr] a dri b , (2.3) 



where again h° b is the metric on the unit hyperboloid. In four dimensions there are two additional conditions 
that must be imposed [3, 26, 27]. First, we require h^ = —2ah^. This rules out the possibility of the so-called 
"spi super-translations" [3]. Second, the function a must be even under the natural inversion mapping induced 
by the Minkowskian transformation x a — ► —x a , ruling out the so-called "logarithmic super-translations." These 
conditions ensure that the asymptotic symmetry group at spatial infinity is precisely the Poincare group. See 
[3, 4, 5, 25, 28, 27] for further details. 

Given the form of the metric (2.2), it is natural to decompose the Einstein equations using the (outward- 
pointing) unit vector n a = N^ 1 5 a normal to the constant p surface TC P with induced metric h ab = g ab — n a n b . 
Denoting projection onto TC P by _L, this yields 

±(R ab ) = K ab + V a a b -a a a b - £ n K ab - KK ab + 2K a c K db , (2.4) 

±(R ac n c ) = V h K ab -V a K = -V\ ab , (2.5) 

R ab n a n b = -£ n K-K ab K ab +(v b a b -a b a b ^j, (2.6) 

where lZ ab is the intrinsic Ricci tensor on Ti p , T> a is the (torsion- free) covariant derivative compatible with the 
induced metric, and £ v denotes the Lie derivative along a vector field v a . The 'acceleration' a b and extrinsic 
curvature K ab are 

a b = £ n n b K ab = l -£ n h ab . (2.7) 

One may take the fundamental variables to be the metric h ab on TC P and its conjugate momentum ir ab = 
h ab K -K ab . 

We now discuss the asymptotic expansions of equations (2.4) - (2.6), taking care to distinguish between the 
d = 4 and d > 4 cases. 

Asymptotic Expansion: d = 4 

The leading term in the asymptotic expansion of the field equation (2.4) implies that the three dimensional 



metric h° b is a solution of 



^ = 2«. (2.< 



By inserting the expansion (2.3) into the field equations, Beig [26] showed that the first order Einstein equations 
are identically satisfied if 

D 2 a + 3a = 0. (2.9) 

Here we have introduced the (torsion-free) covariant derivative D a on the hyperboloid compatible with the 
metric h°^ . Beig [26] also showed that the second order equations may be written in the form 

h {2)a = X2 a 2 ^a a a a , (2.10) 

D b h^ b = 16aa a + 2a b a b a , (2.11) 

D 2 h { V - 2 h { 2 = 6 a c a c h ( 2 + 8 a a a b + 14 aa ab - 18 a 2 h { ^ + 2 a ac a c b + 2 a abc a c , (2.12) 

where a a = D a a, a ab = D b D a a, and a abc = D c D b D a a . For further details the reader is referred to [25, 26, 23]. 

Asymptotic Expansion: d > 4 

It is straightforward to carry out calculations analogous to those of [26] for d > 4. The details are presented in 
appendix C.l; a discussion can also be found in [17]. The equations obtained from the asymptotic expansion of 
(2.4) are 

<? = {d-2)h® (2.13) 



*2> = DA'-id-l^S-^K^ + ^hU (2.14) 

*S> = <.d-2)(h>$-\h™h l 2) . (2-15) 



2 
and the expansion of (2.5) yields 



EPhSH-DahP = 2^-|Wa (2.16) 

D b h%-D a hW = 0. (2.17) 



The function a satisfies an equation similar to (2.9), which is obtained by taking the trace of (2.14) with hfy 
and using the covariant divergence of (2.16) to simplify the result. This gives 

D 2 a + {d _ 3)id _ 1)a+ (d- 3 )(d- 4 ) h (i) = (2 . 18) 

In four dimensions the equations for the second order term h^ b receive contributions that are quadratic in hV 
and a, because these terms all appear at the same order in the asymptotic expansions. This is no longer the 
case in d > 4, and as a result the analog of (2.10) is 

h {2) = 0. (2.19) 

Note that in the asymptotic expansion of Minkowski space the metric h b is a metric for the maximally 
symmetric hyperboloid, i.e., for de Sitter space. We will only require h 1 to be a general Einstein metric 
satisfying (2.13). However, we will continue to refer to 7i as the 'hyperboloid', for convenience. 

2.2 Electric Part of the Weyl Tensor 

The electric part of the Weyl tensor figures prominently in several of our results. It is defined as 

Eab ■= C acbd n c n d . (2.20) 

On a constant p surface 7i p this can be given a purely geometric expression in terms of the intrinsic and extrinsic 
curvatures, the acceleration a b , and the covariant derivative T> a 

Eab = (ir^O (i~[ hab£nK ~ £nKab ) + (^\l KaCRbc (2 - 21) 

+— !— h ab K cd K cd + — L- K (K ab - -J— h ab k] - — L- (n ab - -i- h ab K 



d-1"" d-2 ~ y~™ d-l uu J d-2 V d-\ 

fd-3\ / 1 , m r 1 , 

+ I -TTo I \ ~ J^l -a a a b + — — h ab a c a 

It is straightforward to verify that E ab is traceless and divergenceless. In d > 4 dimensions the asymptotic 
expansion takes the form 

E ab = p 3 - d E$ + p 2 - d E® + Oip 1 -') . (2.22) 

However, there are differences between the terms in this expansion in the d = 4 and d > 4 cases. In four 
dimensions they are given by 

ES = -<r*-h<$a (2.23) 

Eab = -h a 2 b ] +aa ab + bh^a 2 -2a a a b + h^a c a c , (2.24) 



while in d > 4 dimensions the terms in the expansion are 

*ff = - (< -^-\ w -.,-M-iliff. (2-25) 

£f , = Jrf-3^^. (226) 

The structure of these terms and the role they play in the construction of conserved charges are discussed in 
detail in section 4. 

In four dimensions the second order term (2.24) satisfies a non-trivial identity 

*% ~ *&2 = £c d (aD c P b <], (2.27) 

where (5 ab is the second order magnetic part of the Weyl tensor [26, 3, 4, 23]. This identity is useful in establishing 
the equivalence between our expressions for the Lorentz charges and the results of Ashtekar and Hanson [3]. 

3 A Variational Principle and the Boundary Stress Tensor 

It was shown in [10] that a good variational principle for gravity with asymptotically flat boundary conditions 
in d > 4 dimensions is given by the action 

where K := h ab K ab and K ab is defined to satisfy 

n ab = k ab k - k a c k cb . (3.2) 

This equation (3.2) admits more than one solution for K ab - we choose the solution that asymptotes to the 
extrinsic curvature of the boundary of Minkowski space as dAi is taken to infinity. As described in [10], 
the boundary term can be motivated from the heuristic idea that one should subtract off a "background" 
divergence. The defining equation (3.2), arrived at via the Gauss-Codazzi equations, identifies K ab as the part 
of the extrinsic curvature that is fixed by the kinematics of the theory. A similar boundary term has been used 
to construct variational principles for asymptotically linear dilaton spacetimes [29, 30] and asymptotically plane 
wave spacetimes [31]. 

Varying the action and imposing the equations of motion yields 

ss = ttt; I " / ~ K (** - * ab + Aab ) Sh °» ( 3 - 3 ) 

16?rG J dM V ) 

where 7r ab := h ab K — K ab , in analogy with 7r at> , and A ab represents a number of 'extra' terms. An expression for 
A ab is given in appendix A. The boundary stress tensor is defined as the functional derivative of the on-shell 
action with respect to h ab 

T ab ._ _? == ML = J_ f^ab _ ^ab + A ab\ (3^ 

yj—h oh ab 87rG V / 

The first two terms in (3.4) were considered in [10, 23] 

Tab" = ^q (Kab ~ K ab ) . (3.5) 

These terms are conserved in the sense that T> b T^ = 0, and the generators of asymptotic symmetries con- 
structed using T^ w agree with the Ashtekar-Hansen charges. However, the last term in (3.4) was overlooked in 
previous investigations of the stress tensor. In the next section we show that this does not alter the conclusions 
of [10, 23]. 

6 



4 Conserved Charges 

By the general arguments given in [15, 10, 22], the boundary stress tensor leads to conserved charges that 
generate the expected asymptotic symmetries. The generator of diffeomorphisms along an asymptotic Killing 
field £ a is 



Q[£] := lim / J-h Cp u a c T ab e, (4.1) 

where each C p is a Cauchy surface within a constant p hypersurface TL P such that C = linip^oo C p is a Cauchy 
surface in the boundary hyperboloid TC, and Uq is the unit normal to C p in TC p . Explicit computations in 
[10, 23] showed that for d = 4 the quantity 



Q[£] := lim ^/-h Cp u a c T^e, (4.2) 

computed without the term A a }>, agrees precisely with the Ashtekar-Hansen conserved quantities [3]. In section 
4.1 we show that the extra terms in T ab do not change this result. The conserved charges do not receive 
contributions from A ab , despite the fact that in d = 4 these terms are non-zero at the relevant order in the 
asymptotic expansion. The proof relies on a detailed calculation of A ab in appendix B. The d > 4 case is 
treated in section 4.2. A calculation of A ab in appendix C shows that the extra terms are identically zero for 
d > 4, up to terms that vanish too fast in the limit p — > oo to contribute to conserved charges. 

Asymptotic Killing Fields 

Recall that Minkowski space Killing fields naturally divide into two classes: translations and Lorentz transfor- 
mations. In an asymptotically flat spacetime there is a one-to-one correspondence between these classes and 
the first two terms in the asymptotic expansion of a general asymptotic Killing field £ a . In dimension d > 4 
this expansion takes the form 

e = itv) + p- l e { i) + o{ P - 2 ). (4.3) 

Spatial rotations and boosts are associated with the leading order term in the expansion. They correspond to 
Killing vectors £" 0) of TC and therefore satisfy Dr a £(%\ = 0. Translations appear at next-to-leading order in 
(4.3), and are associated with conformal Killing vectors £?U that satisfy £>( a £ (1) M = f(r])h ab for some non-zero 
function / on TC. Terms of order p~ 2 and higher are not relevant for our purposes - the asymptotic expansion 
of the stress tensor implies that they do not contribute to the conserved charges. 

In d = 4 the surface TC at spatial infinity is the unit hyperboloid, so the meaning of asymptotic translations 
and Lorentz transformations is clear. However, in dimension d > 4 we only require hS to be a general Einstein 
metric satisfying (2.13). In that case the properties described above, which depend on the isometries and 
conformal isometries of the metric at spatial infinity, define what we mean by "translations" and "Lorentz 
transformations." 

Potentials for the Stress Tensor 

We now pause to establish a few useful results concerning potentials for transverse, symmetric tensors. The 
existence of these potentials plays an important role in the construction of the conserved charges. Much of our 
treatment follows that of [3, 4, 26]. 

A tensor Q ab is said to admit a scalar potential a if 

9 ab = D a D b a - h$ D 2 a - (d - 2) a h$ , (4.4) 



where the metric h ab satisfies (2.13). Taking the divergence of 9 ab and commuting covariant derivatives shows 

;" 0) is a Killing vector of h° b then the current 6 ab ^ ()) 



that it is conserved: D a 6 ab = 0. If £? 0) is a Killing vector of h°, then the current 6 ab $ 0) can be expressed as 
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the divergence of an anti-symmetric tensor 

9 ab e m =D b (2^ 0) [ b D a] a + aD [b C%) . (4.5) 

Now consider a Cauchy surface C GTC, with timelike unit vector u? Q y The component of the current along u? -> 
becomes a total derivative on C 

«(O)0«6£(O) = u a {0) D b (2ti%D a] a + aD [b Z%) (4.6) 

= K 



2^°\ a D b] a + aD [a C%)u b {0) 



This vanishes when integrated over C, so that currents of the form 6 ab £?„.. do not contribute to the conserved 
charge associated with £" 0) . 

Similarly, a tensor t ab is said to admit a symmetric, transverse tensor potential if 

t ab = D 2 lab + 2U^ bdl cd (4.7) 

for some symmetric 7 a f, with D a 'j ab = 0. The tensor t ab is conserved, and for £" 0) a Killing field of h^ the 
current t a fe £, b 0) is the divergence of an anti-symmetric tensor 

t ab ef 0) = 2 D a (e {0) D [alb]c + ~, c[a D b] e m ) . (4.8) 

A current of this form does not contribute to conserved charges, because the component along u? -, is a total 
divergence on C. The proofs of these statements, which are more involved than in the case of the scalar 
potential, allow for metrics with non-vanishing Weyl curvature. 

4.1 d = 4 

In this section we show that A ab does not contribute to the conserved charges in four dimensions. That is, 

Am = ^ ^G X/^A* ? = (4.9) 

for any asymptotic Killing field £ a . We show that AQ[£] is zero using the asymptotic expansion of A ab . The 
main result of this calculation, carried out in detail in appendix B, is equation (B.45) 

A ab =: pAg + A«+p- 1 A2>+0(,- a ) 



1 



9 o- c a c h^ - 29 a a a b + 63 oo ab + 24 a ap a b p - 5 a cd a cd f$ + 45 a 2 h$ 



3 o- mnp a mnp h ab + 9 o- pq ( a o~ P9 b) ~ 3 cr pq a pq ( ab ) - 2 a e a e(yab ) 



+ 0(p~ 2 ). (4.10) 



Three remarks are in order here. First, both A^° b and A),V are identically zero. Second, all of the terms in A^ b 
involving h^ b cancel, so that it depends only on a and its derivatives. Third, A^ is traceless: h^°' ab A^ b = 0. 

In [23] it was shown that A7r a f, := Tt ab — Ti ab is conserved. From general arguments given in [10, 22] we expect 
the full T ab to be conserved as well. It follows that A ab should be independently conserved. Using various 

(2) 

identities from the appendix D one can verify that the divergence of A ab is indeed zero: D a A ab = D a A ab = 0. 
Thus, any contributions to the charges from (4.9) would also be conserved. Since A ab = there are no 
contributions to the charges associated with asymptotic translations. There may, however, be a contribution to 
the Lorentz charges from the non- vanishing term A^ 



c 



AQ[£] = ^ / y/-t® «8« A i? & - (4- 11 ) 



SirG r y c (0) ab ' (0) ' 



where the integral is performed over a cut C of the unit hyperboloid 7i. In the rest of this section we will only 
be concerned with Lorentz charges constructed from asymptotic Killing fields with £ { a 0) 7^ 0. 

(2) 
Most of the terms in A ab can be expressed in terms of the scalar and tensor potentials described above. An 

explicit calculation shows that 

4 A^ } = 6 ab + t ab + - lab + — j ab - 4 k ab (4.12) 

where 6 ab and t ab were defined in (4.4) and (4.7), respectively, and 

a = --a mn a -a a m , (4.13) 



1 
2 ~ ab -"■'<" 2 



lah = 6(2 aa ab + ~ a 2 h { 2 + a a c a cb - - a mn a mn h$ ) , (4.14) 



2 ...ff+^ + i^"" 



fe a6 = -TT^fcC+^^+o^C- ( 4 -!5) 



j afe - 2( ( j ffa6 + 2a 2 ftS? + ^6 - cr k a k hfA . (4.16) 



Note that k ab , j ab , and j ab are all divergence free. Furthermore, ^ ab and j^ themselves admit similar potentials 

[23] 

\labi\o) = D a (il 0) D [a k b]c + k c[a D b] e {0) )+jabe { oy (4.17) 

jabtfo) = D a (a 2 D [a ^%-4aa [a ^ . (4.18) 



'(0) 

Thus, with the possible exception of k ab , the terms in (4.12) cannot contribute to the conserved charges. 

We will now show that the k ab term does not contribute to the conserved charges. First, consider the case 
of spatial rotations. It is convenient to introduce coordinates (r, 0, (ft) on the unit hyperboloid H. so that the 
metric takes the standard form 

h i a ° b ) dx a dx b = -dr 2 + cosh 2 r {d8 2 + sin 2 6 deft 2 ) . (4.19) 

Without loss of generality the cut C is taken to be the surface r = 0, or "neck", of the hyperboloid, which is 
just the unit two-sphere. In that case the timelike vector u^ 0) points in the r-direction, and for any rotational 
Killing field £" 0) we have 

<o) *2? & = • ( 4 - 2 °) 



(0) n ab MO) 

. 4" Oriri tti 1/1 1 ) 1 i m 1 ran rv\ r fho l r~i -^/nrrvo I i^\t .t .t C 

(0) 



Thus, the contribution to the conserved charge from the k ab term in (4.12) is given by the integral of cr T a" a £" 



over the neck of the hyperboloid 

Qk = j \a T a a ^ 0) . (4.21) 

Since a has even parity on Ti it is convenient to decompose it as 

a = T + Y + + T_Y_ (4.22) 

where Y± schematically denotes even (odd) parity spherical harmonics on S 2 , and T± is an even (odd) parity 
function of the coordinate t alone. Note that T_ is identically zero at r = 0. Using this decomposition we have 

a T = f + Y + + TLF_ (4.23) 

and 

a a ^ 0) =T + Y +i + T_Y^, (4.24) 
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where Y±£ = £? 0) D a Y±. Thus 



Qk = (T_T + ) / Y +i Y_ 

! S 2 



(4.25) 



But Y + £ has even parity when £2L is a spatial rotation, so the integrand in (4.25) is odd and Qk vanishes. 

A similar argument holds for boosts. On the r = slice of the hyperboloid, a boost Killing field £ { a 0) can be 
expressed as 2 



£?o)=0«(o) 



where g is a function on5 2 with odd parity [32]. Thus, 



s* 



Qk ~- I u (ff) kab S(Q) — / u (Q)kab U (Q) 9 — / ^tt 9 • 



S-' 



s2 



(4.26) 



(4.27) 



From the definition (4.15) we see that k TT is a sum of squares, and therefore has even parity on S 2 . But since 
g has odd parity on S 2 , the integrand in (4.27) is again odd and Qk vanishes. This completes the proof that 
the extra term A ab in the stress tensor does not contribute to the conserved charges in four dimensions. 

Asymptotic Expansion of the Full Stress Tensor 

We close this section by summarizing our results for the boundary stress tensor in four dimensions. The 
asymptotic expansion of the full stress tensor is 



Tab 



1 



8ttG 



Air 



(i) 

ab 



A«S?+A» 



+ 



(4.28) 



The large p expansion of A7r a fe was performed in detail in appendix B of [23], where it was shown that 



Air 



An 



(i) 

ab 

(2) 
ab 



(0) 



Cab + 0- h ab 



-E 



(1) 

ab 



M 



a 2 + a c a c + ^ a cd a cd ) ti% + ti£ + 2 a a o b + aa ab + a a c a cb 



,(2) 



-E a 2 b ] + \lab 



(4.29) 



(4.30) 



This expansion can be used to extract expressions for the generators of asymptotic symmetries from (4.1). The 
term A ab does not contribute, so the relevant terms are 



Qlt] 



l 



lim 
8irG p-^-oojQ 



d x 



-h 



(0) 

c 



.-(I) tb 



p(1) tb 



E 



(2) 
ab 



aE 



(l) 

ab 



6 lab ) 3°) 



(4.31) 



P U %) E ab £(0) + U (0) ^afe £(1) + U (0) 

Two of the terms in this expression vanish when integrated over C. The integral of the first term vanishes 
because £|m is a Killing vector on the hyperboloid and E^ b takes the form (4.4), with scalar potential a. The 
integral of the term involving ^ a b also vanishes, by virtue of (4.17). The remaining terms are 



i 



8ttG 



d x 



c 



-h 



(0) 

c 



\0) ^ab ?(1) + U l( 



(0) \^ab 



(1) 



°e Me 



'(0) 



(4.32) 



where the <?E^ b term is due to sub- leading terms in the asymptotic expansions of \J —he and u a . The charges 
for asymptotic translations come from the first term in the integrand, and the Lorentz charges come from the 
second term. The equivalence between the stress tensor and Ashtekar-Hanson definition of the Lorentz charges 
follows from the identity (2.27), applied to the second term in (4.32) [23]. 



2 On a different cut the vector will have other components. 
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4.2 d>A 

The asymptotic expansion of A a b in dimension d > 4 takes the form 

A ab = P A<°> + p~^ A« + p"^ 3 ) A« + 0(p-^) . (4.33) 

In appendix C we show that A^° 6 = A^ = A^ 2 fe = 0, so the stress tensor takes the form 

Tab = ^—^ i^ab - ^ab) 

This is a pleasing result because it is proportional to E ao , without the additional terms that appeared in four 
dimensions. The relation between the stress tensor and the electric part of the Weyl tensor is 

T - = -87G^3 £ - + 0( ^»- (4 ' 35) 

The expression (4.1) for the conserved charges becomes 

m = ~ sTg (a=3) A™ \f~ 2x yf~^ [ pu h E 2 $>) + S) ^ & + u (o ^? *)] • ( 4 - 36 ) 

The first term is proportional to p, and must vanish for the p — > oo limit to exist. In d = 4 the analogous term 
vanishes because E^ admits a simple scalar potential. The proof is more complicated for d > 4, but proceeds 
in essentially the same way. The term E]J can be expressed in terms of one tensor and two scalar potentials as 

,(i) - r , d-4 /. „ „ fd-l 



Kb = - e ab[<r] ~ — K& " <W>» (1) ] - 2 I — - 1 a6 [a]J , (4.37) 

with the potential for t a b given by 

t^c-^ 1 ' -2(^)4?- (4.38) 

It follows that the flux across C of the current E^ ^ 0) vanishes. Taking the p — » oo limit gives the final 
expression for the conserved charges in dimension greater than four 



1 1 



8ttG (d-3) 



r* 



2„J JO) 



GE] = -^ t^t ^* V ~ h c <o) Kl' &) + <o) Kt' e o) • (4-39) 



(!) t& i „.a F (2) t 6 



We remind the reader that the above expressions hold when h°£ is any Einstein metric satisfying (2.13). As in 
the case d = 4, the charges for asymptotic translations come from the first term and the Lorentz charges come 
from the second term. However, for d > 4 the sub-leading terms in the expansions of \/—hc p and Uq do not 
contribute to (4.39) as they begin at too high of an order in p _1 . 

5 Examples 

It is useful to illustrate the results of the previous sections with a few concrete examples. In the first example 
we determine the stress tensor and conserved charges for the Kerr solution in four dimensions. Then, in the 
second example, we explicitly construct four dimensional spacetimes with AL 7^ 0. 
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5.1 The Kerr Black Hole 

The line element for the Kerr solution in four dimensions can be written as 

ds 2 = ^-dr 2 + Zd9 2 + (r 2 + a 2 ) sin 2 9 d4> 2 - dt 2 + 2M ^ r [ a sin 2 9 d<j) - dtf . (5.1) 

The parameters o and M are related to the mass and angular momentum, and the functions S and A are given 
by 

E = r 2 + a 2 cos 2 {8) A = r 2 - 2M Gr + a 2 . (5.2) 

The line element must be expressed in the Beig-Schmidt form (2.2) before the results of the previous sections can 
be applied. This is accomplished by a series of coordinate redefinitions, as described in [25]. In the Beig-Schmidt 
coordinates, h°^ is the standard metric on the hyperboloid 

h { ^dx a dx h = -dr 2 + cosh 2 r (d9 2 + sin 2 9 dcj) 2 ) (5.3) 

and the function a is given by 

a = M G cosh 2t sech r . (5.4) 

Notice that this is the only admissible solution of (2.9) that depends only on r 3 . The first-order terms in the 
expansion of the metric are h a l = —2ahJ ) , as required by the asymptotically flat boundary conditions. The 
second order terms h^ are a bit more complicated. Their particular form is not too illuminating; it is sufficient 
to note that h~l is the first term in the asymptotic expansion that depends on the rotation parameter a. 

The Kerr solution has two non-zero conserved charges: the mass associated with translations dt and the 
angular momentum associated with rotations d^. Calculating the charges from (4.32) requires the t-t and T-(f> 
components of the electric part of the Weyl tensor. Specifically, the terms we need are 

Eg = -2MG sech 3 r E^) =0 E%} = 3 a M G sin 2 9 sech 2 r . (5.5) 

The mass is obtained by evaluating the first term in (4.32) for the vector field 

f° = p' 1 cosh(r) 5 a T + . . . . (5.6) 

The cut C is taken to be the r = surface of the unit hyperboloid, with timelike normal u? \ = 5^ T . Then the 
conserved charge is 

Q dT = -^ />x V^ sin 0u[ o) 4^(1) 

1 /"27T fK 

= — — dcp d9 cosh 2 r sine • 2MG sech 3 r • cosh r (5.7) 

8vrG J J 

= M . 
The vector field for rotations is just £ a = 5 a i . Using this in the second term of (4.32) gives 

1 



q* = -^g i c Sx v - h c -h [m - **%) % 



i 



8ttG 
-aM 



/"27T /*7T 

/ d<j> / d9 cosh 2 Tsin9 -3a MGsech 2 T sin 2 9 (5.8) 

Jo Jo 



These are the expected results for the mass and angular momentum of the Kerr solution. 

v^ for the Kerr solution. Using the exp± c ,, k,:,,, ,*,, ,, iih 



It is interesting to calculate the extra terms A' 2 , for the Kerr solution. Using the expressions for h°l and 



a, all of the components are found to vanish 



Ai? = . (5.9) 



ah 

This result also applies for the Schwarzschild and Reissner- Nordstrom solutions. In these cases the results of 
section 4.1 are superfluous, since the boundary stress tensor reduces to the form considered in [23]. 



3 A second solution is not even under r — > — r, as required by the boundary conditions. 
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5.2 Spacetimes with A$ ^ 

In the previous section we saw that the extra terms vanish for a number of familiar four dimensional spacetimes. 
However, it is easy to demonstrate the existence of asymptotically flat spacetimes for which the extra terms 
do not vanish. The condition (2.9) is a linear equation, so new solutions may be obtained by superposition. 
The result (4.12), on the other hand, is non-linear, which means that these new solutions will generally have 
non- vanishing A^ b . 

An explicit example is a spacetime consisting of two black holes related by a boost. For simplicity we will 
work with the Schwarzschild solution. Since this is recovered from the a — > limit of (5.1) it follows that a 
for the Schwarzschild metric is also given by (5.4). First, consider a single black hole that is boosted along the 
z-axis by an amount (3. The function a becomes 

a b = - q MG (5.10) 

with q given by 

q = cosh/3 sinhr + sinh/3 coshr cos9 . (5-11) 

The boost simply maps one solution of (2.9) to another, and does not change the fact that the extra terms 
vanish. Now consider a spacetime containing two black holes related by the boost described above. The full 
metric for this spacetime may be complicated. But the first order term in the asymptotic expansion satisfies a 
linear equation, so it is given by the sum of (5.4) and (5.10) 

a = a + a b . (5.12) 

As expected, the extra terms do not vanish for this solution. Assuming that the boost parameter is small, the 
first non-zero terms appear at second order in (5 

A&> = P 2 F 1 (r,6) (5.13) 

A$ = (3 2 F 2 (t,9) sin20 tanhr (5.14) 

Ag = (3 2 F 3 (r,9) (5.15) 

Ag> = P 2 F A (r,9) (5.16) 

4? = A S = (5.17) 

where the F n (r,9) are polynomials in cos 28 and cosh2r. The result is traceless and diver genceless, and the 9 
and r dependence of the first two terms is such that they do not contribute to the Lorentz charges. 

The functions (5.10) and (5.12) contain all the information needed to determine the energy and momentum 
of the solution described above. The energy and momentum for the Schwarzschild solution are the same as for 
Kerr, and for the boosted black hole we have 

E = M cosh/? p z = M sinh/3 p x = p y = . (5.18) 

These charges are obtained from the first term in (4.32), which is linear in a. Therefore the energy and 
momentum for solutions like (5.12) are given by the sum of the energy and momentum of the constituent 
solutions. 

6 Conclusion 

The boundary stress tensor (3.4) obtained from the action (3.1) can be used to define a conserved charge for 
any asymptotic Killing field £ a 

Q[£] = g^ fd d - 2 xy^h^u a T ab e , (6.1) 
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where (C, he) is a Cauchy surface in (Ti., h), and u a is a timelike unit vector in TC normal to C. Such charges 
necessarily generate the asymptotic symmetries of the theory. 

In [23] it was shown a slightly different expression, in which T ab in (6.1) is replaced by 

Tfr = ^ (tab ~ * ab ) , (6.2) 

agrees with other definitions of conserved charges in the literature. However, the difference between A ao = 
T ao ~ Tit* consists of a number of 'extra' terms whose asymptotic expansion takes the form 

A ab = p*- d A« + p 3 - rf A$ + 0(p"- d ) . (6.3) 

The first two terms in this expansion appear at just the right orders to contribute to the conserved charges 
associated with asymptotic translations and Lorentz transformations, respectively. This observation creates a 
certain tension with the results of [10, 23], which did not take these extra terms into account. 

Our main result was to resolve this issue by proving that the extra terms do not contribute to the conserved 
charges. In dimension d > 4 the first two terms in (6.3) simply vanish, as a result of a remarkable series of 
cancelations. The proof in four dimensions is more subtle. In that case the first order term vanishes, but the 
second order term may be non-zero. However, in section 4.1 we showed that A a 2 b takes the form (4.12), which 
implies that it cannot contribute to the conserved charges. The remaining terms in T ao are of the form (6.2), 
so that the analysis of [10, 23] applies without modification. We find that (6.1) can be expressed in terms of 
the electric part of the Weyl tensor. In four dimensions the charges are given by 



Q® - "i^5 If^-hW ufo £$ £(i) + <0) OS? " <^£>) 4,) , (6.4) 



1 h^J /.'"■' 
<kG Jc 

while in dimension greater than four they are 

1 1 



8vrG (d-3) 



c 



Q® - -^ TT~^ I d^x J-hf \uU E® & + u? 0) E&> fa . (6.5) 



The explicit comparisons made in [23] show that these charges are equivalent to both the ADM construction 
[1] and (in d = 4) the covariant charges defined by Ashtekar and Hanson [3]. 
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Appendices 

Throughout the body of the paper indices a,b,c from the beginning of the alphabet are used for all tensors. 
This applies to d dimensional spacetime indices, d — 1 dimensional indices on 7i, and d — 1 dimensional indices 
on a Cauchy surface C C TC. In the following sections tensor indices range into the middle and upper parts of 
the alphabet. This is due to the large number of indexed quantities appearing in some of the calculations. All 
indices should be interpreted in the same manner as the indices used in the main text. 
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A Explicit Form of A 



ab 



In this appendix we derive an expression for the extra terms A ab that appear in the stress tensor. Varying the 
action (3.1) and imposing the equations of motion yields 



5S 



1 



-h 



^ab _ ^b^ + 2 £ab} g hab _ 2h ab 5K ab 



!6vrG j dM 
To simplify this we must compute 5k ab . From the defining equation for K (3.2) we have 



where following [10] we define 



L ab cd 5k cd = dn ab + (k ab k cd - k a c k b d ) 5h cd 



cd ucd i 



L ab ca = h ca K ab + 5 c a 5 d b k - 5 c a K d - 5 c b K d . 
Inserting equation (A. 2) into (A.l) we obtain 
1 



5S 



-h 



^ab ~ab , {sab njcd 



ab if a if b 



^ao _ % ao + Ra0 _ 2L ca KcdR ao _ K c a K d ° 5Il ab - 2L ab 51Z, 



16?rG j dM 
where we define -k = h K — k and where 

lab _ h mn /^-lx ab 
\ 'van ' 

with (i _1 ) ab defined to be the inverse of the tensor on the left hand side of (A. 2), i.e., 

5k ttb = {L- 1 ) ab * 

We can use the fact that 



sua - ( K k ij k cd - ki c kj d ) si,,,. 



-i 3 \^^ij j 



8Ki 



1 



.kU 



1 



Ml 



kU 



brViVjShu - -h^VkVrfhij + h Kl V k V {i 5h j)h 



(A.l) 

(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 
(A.7) 



where T>i is the torsion- free covariant derivative compatible with h ab . We now perform integrations by parts to 
write 



SS 



1 



WttG 



c)M 



n ab _ ^ab + jsab _ 2 jcd ( K cd K ab - K c a K d b 



+V 2 L ab + h ab V k Vil kl - V k (v a L kb + V b L ka 



5h ab . 



(A.l 



The boundary stress-tensor for pure gravity can be readily obtained from this equation 

rpab 

where A ab is defined to be 



2 OO 1 I ab ~ab , A ab 

IT — 7T + L\ 



-h 5h ab 8ttG 



(A. 



A a6 = K ab - 2L cd (k cd K ab - K ca k bd ) + V 2 L ab + KbVkViL* 1 - 2V k V {a L k b) . (A.10) 

Some explicit results were given in [33]. The result (A. 10) is analyzed in detail for the d = 4 and d > 4 cases 
in sections B and C, respectively. 
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B Stress Tensor Calculation d — 4 

This appendix contains details of the boundary stress calculations in d = 4. In particular, we calculate A a ^ as 
an expansion in inverse powers of p. We work with zero NUT charge (i.e., we use h ab = —2ah ab ) and make use 
of asymptotic field equations (2.9-2.12) as needed. The calculation is organized as follows: in section B.l we 
calculate the inverse of L. We then perform p~ l expansions of L and L M^d in section B.2. In section B.3 
we calculate p~ l expansions of the derivative terms needed in (A. 10). In section B.4 we collect various results 
and present a final expression for A ab . Equation (B.45) is the main result of this appendix. Throughout we 
make use of identities from appendix D. 



B.l The Inverse of L ab cd 

Let us begin by calculating the inverse of L a f, . The form of the expansions for L and its inverse are 



L (0) + I L (1) + 1 L (2) + 



(L- 1 : 



(L 



P 

-n(o) 



P' 



Av 



-i^(i) 



+ 



(L- 



-1^(2) 



+ ... 



(B.l) 
(B.2) 



with the coefficients in the expansion of L 1 given by 

~1\(°) ,kl r(l) mn f r -l\(°) 



(L 



l^( 1 ) pq 

ij 



p)% H zw H ™ir 



i.i 



L- 



1\(0) kl t- (2) mn (t-1\(°) 



kl 



L~ 



pq 



pq 



(L-^i^LUki™ 



(L- 1 ) 



1\(0) 



pq 



(B.3) 

(B.4) 



The part of L 1 relevant for our calculations is actually W 3 '{L l )ij ■ As a first step we calculate hP> l:> (L x ) ij pq . 
Using the expression for (L _1 V 0) (equation (B.7) of [10]) and equation (B.3) we obtain 

h (0) ij ( L ~l)(2) ..pq = _P_ h (0) M L (2) m n ^-1^(0) ^pq 



4 

P 

"4 



+ '4 h<® ab L^ ab cd (L-') {0) J l L^ M mn (L- 1 )^ mn Pq 



(B.5) 



The first term in this equation is linear in Z/ 2 ', while the second term is quadratic in L' 1 '. After a bit of 
calculation we find 



L 



(2)« 

ij 



P 



W h (w _ /.(D^a) _ {h (m _ h wkp h (i)i )h 



(0) 



+ ^ ) (p(2)_p(l). /l (l)_ /l (2) + /l (l). /,(!)) 

where u^ 1 ' • v^- 1 ' = u^'^v}- . Contracting with h^ ' 13 yields the first term in (B.5) 



(B.6) 



h (P)i3 L &)U 



2 h^ kl p(2) _ 2 pW kl - fc« kl p0-) - /,(2) kl + ^(i) k Ph (i) 



h (0) kl h (l) . ^(1) _ h (0) kl h (2) + ^(0) fe^(l) . ^(1) + 2 ^(1) (l\m\ p (l) k) i 



(B.7) 
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Using h\ ■' = —2ah\, and p\-' = o"y — on\- [23] the expression (B.7) further simplifies to 



h ^)iO L f)mn^ L -l^) mr p q = h {0) Pqf) (2)_ 2p (2) Pq _ h (2) Pq _ 4(7 2 h (0) Pq _ 4aaPq _ 

Now consider the second term in equation (B.5). A short calculation gives 

P_ h (0) a^W^ (L -1 )( 0) Jl L W ki mn (jL -l )( 0) ^ 

- a ab a ab h^ pq + — a 2 h^ pq + 2aa pq + a pl a t q 



Adding the two terms (B.8 and B.9) one arrives at the result 



^'(L- 1 ) 



1\(2) 



. .VI 
*3 
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fe(°)Mp( 2 ) + 2p( 2 )P« + h®™ + — a 2 h^ pq + YLoa™ 



a ab a ab h^ pq + 4a pl a l q 



(Bi 



(B. 



(B.10) 



B.2 Expansions of L ab and L ab M f 



abed 

ab 



In this section we calculate expansions of L and L M abc d- Recall that L and M abc d are defined as 



L 



ab 



,cd/i — 1 \ ab 



M, 



abed 



h c {L~ ) cd 
K a bK c d - K ac K bd . 



(B.ll) 
(B.12) 



Let us first study the expansion of L . In the previous section we calculated various elements needed to write 
such an expansion. Note that 



L 



ab 



h mn {L- 1 ) 



P z 



^(0) mn _ _/j(l) mn }_ (^(2) mn _ ^(1) mp^(l)n\ 

p p 2 \ P ) 



+ 



(L 



-1\(0) ab 



+ p- l {L 



lfr-l\(l) ab 



n ab + p- 2 (L 



2/r-l\(2) ab , 

) mn "!*••■ 



.!^ 



«■() 



+ P 



-2 



<T 



06 



<7 



/,(°)«M+p- 4 



^(0)mn (L -l )( 2) f 



4 ' V 2 

■A (1)mn (L" 1 )«„ at - (h {2)mn -h il)mp h£ )n ) {L- l ) { ^ mn ab 



+ 



(B.13) 



After some calculation we find the following expressions for the second and third elements of the second order 
term of equation (B.13) 



/jfi)™^- 1 )! 1 ), 



ab 



p(a 2 h^ ab + ao ab ) 



h (2)mn_ h (l)m Ph (l)n\ (£-1) 



-1\(0) ab 



1 



.^^^V^x^V 



(B.14) 
(B.15) 



Adding together these expressions with (B.10) we find that (B.13) simplifies to 



L 



ab 



.th® 



ab 



+ P~ 



a' 



ab 



+ <T 



h (0)ab) + P: 



2 p(2)ab_ h (0)ab p (2) _ 3h (2)ab 



+h (0) ab ac(J e + ^2^(0) ab + ^^ab _ ^ CT ^(0) ab + ^ae ff 6 



+ •• 



(B.16) 
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Now we want to eliminate p^ 2 > ab in favor of {h^ 2 ' ab , a} and their derivatives. This is done using equation (B.10) 
of [23] . Applying the equations of motion one can show that 



s(2) _ 1,(2) 



,(0). 



,(°)^2 



^d„ .1,(0) 



P% = h ab ~ "St*** + 2 ^b ~ f%<r A + ™ab + a a c a cb - -o^a^h™. 



Substituting (B.17) in (B.16) and doing some algebra we find 



L 



nil 



.^W 



ab 



+ P~ 



a 



ab 



+ <T 



h (0)ab) + P: 



99 



-a 



2/^(0) ab 



T&cdC 



cdifO) ab 



+6a ac a c b + 18aa ab - h {2) ab + Aa a a b - a c a c h^ 



ab 



+ 



The expansion of M abc( i is relatively straight forward 

Mw - o 2 (h {0) h {0) -h^h {0) ) + o(v {1) h {0) + h {0) 3 (1) -v W h {0) 

Mabcd P K n ab n cd n ac n db > + P\P a b a cd ^ n ab Pcd Pac a bd 



'ab "cd + n a6 Pcd Pac "-fed n ac Pdfe / 
^ \Pab n cd + Pab Pcd + 'l ab P cd P ac n db Pac Pdb n ac Pdb > + • • • • 

Using the expansion (B.18) for L and the expression (B.17) for p ab we find 



fab M _ P h (0) °cd + 5ah^ 

L> Mabcd ~ ^cd 4 



-19aa cd - 10a c p a dp 



1 



+ .. 



n a pq a 



7 ^cd ~ 2c7 ^d - y o- /i cd + 2/i cd 



(B.17) 



(B.18) 



(B.19) 



(B.20) 



The expansions (B.18) and (B.20) are the key results of this section. 

B.3 Calculation of the Derivative Terms 

In this section we calculate the p _1 expansions of the derivative terms. We start by calculating a general 
expression for two derivatives acting on L ab . For the moment we keep all four indices free 



V c V d L mn = dWd®L<® mn + p- 1 (pWDWLW mn + D^D^m " 

+p- 2 (D^Df!® mn + BfD^L^ mn + DfDfm mn 



+D^Dfl^ mn + D^DfL^ 



+ .. 



(B.21) 



In expression (B.21) D*-°> is simply D - the torsion free covariant derivative compatible with the unit hyperboloid 
metric h b . D>- ' and D'- 2 ' respectively denote the first order and the second order corrections to the covariant 
derivative T> away from D. Calculating various terms in equation (B.21) requires some amount of work. Using 
equations (D.18) and (D.19) of appendix D we find the following results: 



D [0] D [0] L (1)'. 



- „mn I _ t (0) mn 

-a dc + Odch K ' 



(B.22) 



2. 



DfDfm- 



l h (0)mn 

2 



(B.23) 
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7. 



D [0) D [0]l(2)mn 

c a 



^a c a d h^™ + ^oo cd h^ 



-\<J pq <? pq dch {0) mn - 9 -a pqc a pq d h^ mn 

+Gv mp dc<T P n + 6v mp d <T P n c + 6cT mp cV P n d + 6v mp a P n dc 

+I8a dc a mn + 18a d a mn c + 18a c a mn d + 18aa mn dc 

-h^ mn . dc 

+4a m dc a n + 4a m d a n c + 4a m c a n d + 4a m a n dc 

-2a pdc a p h^ mn -2a pd a p c h^ m 



(B.24) 



D [0] D W L (l)mn 



a ec a< n h^ )m ^ - a e a e( - n c h { ? )m ^ 



eU cn, d 



_ „mn „ „mn 

a dc a - a d a c 



-2a dc ah^ mn - 2a d a c h^ mn + a^ m c a d ^ + a^ m a d n \ 



(B.25) 



D^D [2] L^ mn 



ih® ' 



-a c a d n^> mn - aa cd 



l(0) mn _|_ _l,(2) mn 
4 



dc 



(B.26) 



6. 



D [l] D Wl(0)mn 



\a^ e h^ mn h^ 



(B.27) 



D^Dfrn 



2a [d a mn c] - hfja e a mn e - 2a e a< n d h^ m \ + 2a c ^ d a^ 



2a e a e h^ mn hfJ 



(B.28) 



Various contractions of the above expressions give the required three final terms in (A. 10). 



Term 1 



h ab V m V n L mn = — [45a 2 - 29a m a™ + 5a mn a mn + 3a mnp a mn ?] h^ } + . 
This result uses the following two identities from the appendix D 



(B.29) 



m P n 

u nm u P 

m P n 
u n u P m 



o^- „mn c\Ji 

2<J mn a — y<7 , 

mnp q m 



(B.30) 
(B.31) 
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Term 2 



V 2 L ab 



§K +'-)+£ 



7T a mnp a p h ab + — a m a /i afe + — a h ab - 8 a a a b 



2 



31 



&mnO' nm h ab + 44 a am a m b + 72 <7cr a fe + 12 a a mn <J brnn + 18 cr e cr a b e 



+ 



The calculation of this term relies on the following identity from appendix D 

_ nem _ vane „ra n , _e_ 1,(0) 

over = <T e (T —a a + a a e nL,L. 



(B.32) 



(B.33) 



Term 3 

The calculation of this term is more involved than the previous two terms. We find 



2£>n£ , (ft£ a ) n 



0) H- 
p 



-3{a ab + ah2>) + 



.2,(0) 



5 51 

t VaOb ~r c^ab ~ ^cr ap a b p + a ap qO q + 3 a p a p h 



(0) 
u6 



,(o) 3 



-21a'h%+4cT mn a™h% "^^) - 4^(.^|6) 
The identities from appendix D that are needed for this calculation are 



+ 



Vp(ab)q& 
a (a Pq(J \pq\b) 



'"' ~ v P q(ab)V P9 ~ 6 crvab ~ hf b > a pq a pq - a ap a p b , 
& P q(aV pq b) ~ 3 <J a o b - cr q (T q (a,b) 

= a pq(aV P9 b) - 2 °aOb ~ ^^abq ~ h ab <J p (J P , 

&(a n b)n = 2 cr ab + 3 ah ab . 



(B.34) 



(B.35) 
(B.36) 

(B.37) 

(B.38) 



B.4 A Final expression for A ab 

In this section we collect the expansions performed in the previous sections and give a final expression for A ab - 
Recall that 

A afe = K ab - 2L cd M abcd + V 2 L ab + KbD^L 1 * 1 - 2V k V {a L k b) . (B.39) 

Various terms in A ab admit the following expansions: 
1. 



K ab = ph nh +a ab -ah nb +- 



Kb ~ h ab a ccr c + 2 a a a b - - h^ b 'a z + aa ab + a a c cr cb 



4 



1 „cd„ h (o) 



+ ... , 



(B.40) 



2. 



- L a0 M 



abed 



Pu(o) 



1 



^: + t^+5^7 -- 



(0) 



1 

Tp 



-a <j pq h {0) 



2^ -^2? +2/$ 



•19ao- crf - 10cr c p a d p 



+ .. 



(B.41) 
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habD k T>iL 



ki 



1 



— [45 a 2 - 2$a m a m + 5a mn a mn + 3a mnp a mnp ] h a ° b > + 



(B.42) 



V 2 L ab = \ (cr/$ + a ab ) + ' 



4p 



o- m npCr ^ /i a6 + — a m a h ab + — a h ab - 8 cr a cr b 



-— o- mn a nm h^' -r 14 (J arn a m b + 72 <7cr ab + 12 <7 a m? V 6mri + 18 a e (rabe 



+ ... , (B.43) 



5. 



2T> n V^ b L a - ) n 



-H<rab + <rh$) +- 



51 

T 



7 Ca^b — — CTO"afe — 11 <?ap&b P + VapqV 9 + 3 <7 p (7 P h ab 



-21 a 2 h^+ 4 a mn a mn h 



(o) _ 3 

ab /i 



u u pq(ab) 



a (a a \pq\b) 



+ ... . 



(B.44) 



Adding these terms we find 
A a6 



; ^ , + °(? 



1 
47 



9 <yc° c h { $ - 29 a a a b + 63 aa afe + 24 a ap a b p - 5 a c ^ cd h$ + 45 a 2 fcj™ 



,(0) 



3a mnp a mnp /i a6 + 9a pq ( a a pq b ) - 3 <J pq <T pq ( a b) ~ 2v e &e(ab) 



+ ... . 



(B.45) 



The expression (B.45) is the key result of this appendix. Notice that all of the h^ b terms have canceled, so that 
A v 1 is determined by a and h 



ab ' 



C Stress Tensor Calculation d > 4 

This appendix collects a number of calculations for asymptotically flat spacetimes in d > 4 dimensions. In 
section C.l we perform the asymptotic expansion of the field equations, then calculate the counterterm K^ 
to the requisite order section C.2. These results are used in sections C.3 and C.4 to obtain an asymptotic 
expression for A a b. Remarkably, A a h vanishes identically up to terms that vanish too rapidly (as p — > oo) to 
contribute to the conserved charges. We close in section C.5 by presenting the complete expression for the 
boundary stress-energy. For convenience we set n = d — 3 in parts of this appendix. 



C.l Asymptotic Expansion: d > 4 

We now perform an asymptotic expansion of the d > 4 field equations to obtain the results (2.13-2.19), analogous 
to those of [26, 25] for d = 4. See also [17] for similar results. 
The expansion of the extrinsic curvature for d > 4 is 



ab 



P~ X K 



h 



(») 



l 



ab — 't a b n d-3 



, to , d— 5 (i) 

°Kb + -^—Kb 



d-A h 



(2) 



2 p 



4-2 



+ o 



n d-l 



(C.l) 
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where for convenience we define P ab := p 1 K ab . The trace of the extrinsic curvature is 



K 



(d-1) If,, » d-3, ( i)\ (d-2),( 2) ,/l 



(C.2) 



where h^ = h^ ab h ( ^ and h^ = h^ ab h { ^. The acceleration is 



a b 



-NV b cr 



0, 



1 



P 



d-3 



D h o 



(C.3) 



h (f>) 

n ab- 



to the appropriate order. As in the main text, D a is the torsion- free covariant derivative on 7i compatible with 
Inserting these expansions into the Einstein equation (2.4) yields 




<>-(<*- 2) « 



+ - 



1 



d-3 



P 



,(D 



(0) 



K%-D a D b a+(d-l)ah% + 



d 



n n ab o n ab 



(C.4) 



n d-2 



n 



(2) 
ab 



{d -2)hS + ^h^ 



+ o 



id-l 



The leading term in this equation contains TZ ab , the Ricci tensor for the metric h b . The sub-leading terms 
depend on the tensors appearing in the asymptotic expansion of TZ ab , defined as 



tt&° = \ {D c D a h^ + D c D b ht ] ~ D c D c h% ] - D a D b h^ 



"ab 9 V a cb 

for m = 1,2. The expansion of equation (2.5) gives 





1 



P 



4-3 



(d - 2) DV - ^— - (D b h^ ab - D a h® 



(C.5) 



(C.6) 



(d - 2) 1 

2 ^2 



Dh (2)ab _ D a h (2) 



(0) 



+ o 



1 



,d-l 



where indices are raised and with respect to the metric h , . The remaining equation (2.6) can be written as 



n - K 2 + K ab K ab = 



(C.7) 



which yields 



= ft(0)_(d_2)(d-l) 

' ^ (1) + 2 (d - 2) (d - 1) a + (d - 2) (d - 4) /i (1) 



+- 



+- 



4-3 



.d-2 



■RW + (d - 2) (d - 3) /i< 2 > 



+ 



vJ-1 



Note that 7£( m ) is defined as the trace of (C.5) 

n (m) =h (0)ab n (rn) _ 

This is noi the same as the term at order m in the asymptotic expansion of 1Z. 
To leading order the Einstein equations become 



(C.8) 



(C.9) 



n\ 



(d - 2) h 



ab 



. 



(CIO) 
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The subleading equations are obtained from (C.4), (C.6) and (C.8) by setting the coefficients of the various 
powers of p" 1 to zero. Note that these equations are not all independent, and their self-consistency imposes 
additional constraints on the asymptotic expansion. For example, taking the trace of (C.4) with respect to h~l 
implies that 



D 2 a + (d-3)(d-l)a + 



(d - 3) (d - 4) 



h<U 



and h {2) 



0, 
0, 



(C.ll) 

(C.12) 



where we have used (C.6) and (C.5). 

The leading term in the Einstein equations has different implications for asymptotically flat boundary 
conditions in d > 4 and d = 4. When d = 4, spatial infinity is a three dimensional surface and so the equation 
w°£ = 2 h° b implies that the metric on TL is locally that of the unit hyperboloid. For d > 4 the metric on TL 
may have a non-vanishing Weyl tensor, so h° b contains degrees of freedom that are not determined by (2.13). 
For simplicity we still refer to TL as the unit hyperboloid, but in fact we only require the condition (2.13). All 
of our results in d > 4 apply for metrics h° b that have arbitrary Weyl curvature. 



C.2 Extrinsic Curvature Counterterm 

We now compute the counter-term Kij to the required order. We employ the asymptotic field equations (C.4- 
C.8), expanded in powers of /? _1 , to simplify results as needed. 
To leading asymptotic order the Einstein equations imply that 



K {0) = ph {0) 



using equations (C.4) and (3.2). Expanding this latter equation we have 

PP ab - h cd P ac P db = K { 2 + 
where P a b = p~ 1 K a b admits the expansion 



, s 7? (1) 7? (2) 

cd£) 6 -7,(0) , l^ab lx -ab ,q 



+ 



P" 



P 



,n+l 



„n+2 



1 



1 



It is straightforward to solve iteratively for P \ . This gives 



ab - n ab ^ n ^ab + n +l ^ab 



(I) 



Pt> + o 



,n+2 



P 



(1) 



(ib 



5(2) 
ab 



1 

n 
1 
n 



K 



K 



(i) 

ab 

(2) 
ab 



1 



2 (n + l) 



nWQ-hV + lhWh® 



ab ab 



ab 



L -^ (2) «-« + ^ (2) « 



2 (n + l) 



(C.13) 



(C.14) 



(C.15) 



(C.16) 

(C.17) 



where indices are raised/lowered and traces taken with respect to h° b . When the Einstein equations are satisfied 
these expressions simplify to 



W = l -D a D b o + \h^ PW = I^ + 1 A (D 



P. 



ab 

(2) 



ah 



,(2) 
'ah 



P® = h^ = 



(C.18) 
(C.19) 



where we recall that h^ 2 > = by the trace (C.12) of the second order term in the Einstein equations. 
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C.3 Computation of the Inverse of L ab cd 
Recall that the inverse of L ab cd is defined by 

O )ij Ljkl = Oj Oj . 

Following [10] we introduce the expansions 

and (L- 1 ) = (L" 1 ) (0) + p^^ {L' 1 ) (1) + p^^ {L' 1 ) (2) 
Using the definition of L _1 we find 

1\(!) pq _ _ (r-l\(°) ..kl (t(1)\ mn( T -l\(0) pq 



and 



L~ 



{L-'f\r 



(l-7 



1.1 



n 



ki 



(L-y 



(l-^^Hl^A ™»a-M (0) «»« 



« 



Defining Qjj = p -^i? we nave from (A. 3) that 



,^^ 



Qi/' = fcHji,. + ,*!" ,)',.' P - <)<" P\ ' - 8TP- 



(k pi) _ Jk pi) 

i 3 3 i 



and upon expansion we obtain 



<i3 

'< 1. 1 



n Wkl 



Pffh ( ° )kl - h%h Wkl + sft] (P W - h w ) - 5f (if - hf l) 



,(* (put) _ , «0 



Q 



ij 



P<?h m " - f$h w » + 8f8 l ] (P (2) - h {2) ) - 8? (jf*> - ft? 1 * 

-5? (Pf )l) - h?» 



(C.20) 

(C.21) 

(C.22) 

(C.23) 

(C.24) 

(C.25) 

(C.26) 

(C.27) 
(C.28) 



From (C.23, C.24) we can compute the expansion of (Q l ) . . kl . However the quantity that we really need is 
Q ab = J 7 h rs Q rs ab = ■\h rs L rs ab = p- x l ah , which upon expansion is 



C°)a6 



Q 



Q Wab 



Q (2)ab 



2(n + l) 
1 


n (n + 1) 
1 



n (n + 1) 



(°)afe 



2 2 2 



P 



2 2 



or, using the Einstein equations (i.e., using (C.18) and (C.19)): 



Q^ ab 



Q Wab 



1 



h 



(°)ab 



Q 



C»o6 



2(n + l) 
1 

2ra(n + l) 

1 
~2(n+l) 



h 



-D a D b a - (n + 1) h (1)(l6 + (n + 2) a^ C ° )a6 + W° )c 
n 2 



(C.29) 
(C.30) 
(C.31) 

(C.32) 
(C.33) 
(C.34) 
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C.4 Finding the Full Variation 

We now need to use (A. 10) to compute A Q & 

A ab = j^ab _ 2 ~ L cd (x cd K ab - K c a K d h \ + V 2 L ab + h ab V k V l L kl - V k (v a L kb + V b L ka ) 

which is more easily done in pieces. 

We begin with the first two terms. We have 

K ab - 2L cd (k cd K ab - K c a K d b ) = i (P ab - 2Q cd (P cd P ab - P c a P d b 



(C.35) 



(C.36) 



The zeroth order term is 

pab _ 2 Qcd (p cd pab _ p^ap^ 

as expected. The first order term is 

pab _ 2 Qcd (p cd pab _ p c apj 



«)) 



S°)ab 



2 h^ d (h^h^ b -6 c %)=0 



2(n + l) 



(i) 



pWab_ 2h M<* '" 



(Vcd 



h%(P (1)ab -2h m « b )+P { cd ) h (0) ™) + 



n (n + 1) 
1 
n (n + 1) 



(n+1) 

Wed 1 h m u m cd _ n + 2 h ^cd + l h m h (o) cd 



2 ( p W *b _ h Wab 



(n+1) 



P' ' ca --p 'h 

2 



tGh W *-6 e *6 d > 



^D«D b * - H + V>«* +{[n + 2)a + ^h W ) h^ b 



n 



The second order term has a similar derivation and so 



pab _ 2 gcd fp cd pab _ p c ap d 



V2) 



1 



n (n + 1) 




(n + 2) ( P (2)ab - h (2)ab 



P (2) - h {2) ) h m < 



(C.37) 



(C.38) 



(C.39) 



where the last equality follows from (C.19). 

We now turn to the last four terms in equation (C.35). Using L ab = pQ ab we obtain for the zeroth order 
term 



\kl 



V z Q a0 + h^VkViQ" 1 - V k P a Q ft0 + V°Q K ' 



\a/\kb i T\b^\ka 



(0) 



1 



2 (n+1) 




Vu (a) ab , uabi 



S°)kl 



D A K"" :U + h^D k Dih mia - D k ( D a h Wkb + D b h ( ° )ka 



(C.40) 



again as expected. 

For the first order term we rewrite equation (C.33) as 



where 



X, 



1 



ab 



n (n + 1) 



gWab = X ab 1 ^ab 

* 2 (n + 1) 

D » D b a _ h m "b + I C( n + 2) a + -//" ) fr v 



/? 



2 2 
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(C.41) 
(C.42) 



Making this separation yields the following advantage 



AW 



in 



D a D b Q Wrs + ' 



D a D b X rs + 



2(n+l) 
1 



h mr D a r^ + ^^j 111 ' 



/;/,■ 



2(n+l) 



(°)fcr n t- (1) s , i< >fcsn t- (1) 



/i ,/ep £>aiV + /i ,KS D a r bk r - D a D b h 



a). 



= AA* rs 

because the second term is identically zero. To see this note that 



m rs , l 



l 



D a D b h" rs + -D a [D r hV s - D s hi 1)r ) + -A, (A/^ - D^r 8 ) - D a D b h Wrs 
= 0. 

There will be three kinds of second-derivative terms that occur. These are 

1 



ab 



D a D b X 



D 2 x ab 



A ( D a X kb 



n (n + 1) 
1 
n 
1 



~ AA ( AL> 6 ct) - \D a D h h il)ah + \d 2 [ (n + 2) a + -// 



n, (i) 
2 



ra (n + 1) 

n 
1 



A ( AAct) - l -D 2 h (1)ab + l -D 2 ((n + 2) a + ^h W ) h Wab 



n (n + 1) 



1 n. nH (1) B , 1 n«- 



n, (i) 



A, A, ( £> a £V 



- A^ a ( D*D b a J - -D k D a ti ' kb + -D a D b i( n + 2)a+-h 



To simplify further we need to commute the derivatives. We have, using (D.2), 

= D a D 2 (D a a) = A [D 2 , A] a + D 2 D 2 a 

= -A (TZ {0)aC c d D d a) + D 2 D 2 a 

= {n + l)D 2 a + D 2 D 2 a, 

= [D 2 , A] D b a + A \D 2 , A] a + D a D b (Act) 

= -n^ e c ca D e D b a - K^ eb c a D c D e o - D c hz^ eb c a D e a \ 

_ D a f n (0)e c Cb De \ + D a D b ^ D 2^ 

= 2 (n + 1) D a D b a - 2R { ° )ebca D c D e <j + D a D b (Act) , 

= [A, D a ] (A Act) + A \D 2 , A] a + A A (Act) 

= -K^ ec c a D e D b (T - K (0)ebca D c D e a - D a (lZ^ ecb D e cj) + D" D 1 ' (Act) 

= 2 (n + 1) A Act + A A (Act) - 1Z ( ° )ebca D c D e a. 



A A Act 



A£> a \D k D b a 



(C.43) 



(C.44) 



(C.45) 



(C.46) 



(C.47) 



(C.48) 



(C.49) 



(C.50) 
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Putting this all together we obtain 

V 2 Q ab + h ab V k V x Q kl - V k (v a Q kb + V b Q ka 
= D 2 X ab + h (0)ab D r D s X rs - V k (v a X kb + V b X 



(i) 



n (n + 1) 



1 ' " -- -,2_ n 2 n 2_M.(°)afe 



/? 



((n + 1)DV + D l D l a)h ' ab + 2(n + l)D a D b a 



-2K l0)ebca D c D e( j + I^D 6 (L>V) 



-- (l (n + 1) D a Z)V + D a D b (D 2 a) - K <0)ebca D c D e a 
. 1 (i)V 1)ffl6 + h i0)ab D r D s h mrs - V k (v a h Wkb + P 6 /^ 



n ^W\ (,(°>a6 na n fe 



+L> 2 (n + 2) a + -K' K' ab - D a D° [(n + 2)a + -h 



n,(i) 



1 



n + 2 



1 



ra (n + 1) 

where (D.2) was used. 

From (C.51) and (C.38) we see that 



D a D b a + -(n + 2)h 
n 2 



Wab 



(n + 2)a + ^h W )h Wab 



A 



ab 



(1) 



K ab _ 2 icd (K cd K ab - K c a K d b ) + V 2 L ab + h ab V k V l L kl 



-V k ( £> a L fc6 + V b L ka 



0. 



(i) 



(C.51) 



(C.52) 



In other words, all terms cancel at first order. This is actually a more general result than we have claimed up to 

n (0) 
l ab ■ 



this point, because no assumptions were made about the the Riemann tensor for the metric h° h . The boundary 



conditions require R ? = (d — 2) h J but leave the other degrees of freedom in the Riemann tensor unspecified. 
Our derivation applies for any metric h°, with the appropriate Ricci tensor, even if it has a non- vanishing Weyl 
curvature. 

For the second order terms, using (C.34) we have 



V a V b Q r 



(2) 



D a D b Q (2)rs + ' 



2(n + l) 



,(°)fcr n r < 2 ) s , h m ks n r (2) ' 

n U a^ bk t n U a V bk 



2(n + l) 

0, 



h D a T bk + h D a T bk - D a D b h 



(C.53) 



where the last line follows because the calculation is identical to that given in eq (C.44) above. Consequently 
(C.39) and (C.53) yield 

0. (C.54) 
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ab 



(2) 



C.5 The d > 4 Boundary Stress Energy 

The boundary stress tensor is given by 

rpab _ 



8vrG 



„a6 ~aft 
7T — 7T 



(C.55) 



27 



up to at least second order in the asymptotic expansion. We now compute the explicit form of T . Using 



n ab _ JT ab = J_ ( h ab (p _ p) _ pab + pab 



we find 



pab 106 p _|_ lad p pab 
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pa 



-,ab 



h a °P + h a °P - P 



ab t 
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(1) 



/j(0)a6 _ ( n + 2 ) h^ ab - h^ ab + (n + 2) h^ ab = 0, 



p{l)ab _ p{l)ab + h (0)ab ( p(l) _ p(l) 



(C.56) 

(C.57) 
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D a D b a - - h (1)c 



+h (0)ab A +2 )a + Vi^. h W + (I D 2 a + l - h^X) 



-a 



h<® 



al> 



n — 1 



h™ 



all 



n 



D a D b a, 



and 



pa 



h abp + h abp _ p, 



ab 



(2) 



p{2)ab _ p(2)ab + ^(0)afe f p{2) _ p(2) 



n+ 1 



/l (2)a6_ 



Putting this together yields 
1 1 



->ab 



(a h^ ab + — h W-b + _J_ D a D b\ 1 

V 2 d-3 / p d+1 



1 d-2 



p d 8?rG 

D A Collection of Useful Identities 
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h (2)ab + 



(C.58) 



(C.59) 



(C.60) 



By commuting derivatives several times, one can prove a number of identities involving a and its derivatives. 
Below we present a collection of such identities in d = 4. All of these identities, in one form or another, are 
used in the main text and in appendices B.2 and B.3. For ease of reference we state our conventions again: 
&abc '■= D c D D a a, (J a bcd '■= D c [D c D D a a, etc. Also, recall that the commutator of two covariant derivatives 
acting on a tensor t a b is 



where 



\D a ,D h ]t c ° 



M 



? (0)e, d 
l abc ie 



Md 

L ab eUc ' 



i&^M-'.M. 



One can show using the equation of motion for a that 



Of 


= 


— 3(7 , 


O ca 


= 


-3cr a 


_ c 
°ac 


= 


—cr a , 


_ rnen 


= 


_ vane _m „n , _e_ (,(0)mn 


pnm 
u pmn u 


= 


pmn r> m 
Vpmn® & Vm" i 


u u apq 


= 


V m Op qa ~ 3 OO a - (T a p a p , 


°pq(aVb) Pq 


= 


Vpqia^b) ~ 3 <7 a 0-£> ~ V Q Vq(ab) , 


mn 
"{a < -'b)mn 


= 


Vmn{aV mn b) ~ 2 C n cr n ( o6 ) ~ 5 0" a <7fc + h^ (T m (J m , 


Gahc 


= 


6 CT Kb + 3(J ab, 



o, 



(a b)n 



2a ab + 3(rh < $, 



(D.l) 
(D.2) 

(D.3) 
(D.4) 
(D.5) 
(D.6) 

(D.7) 
(D.8) 
(D.9) 

(D.10) 

(D.ll) 

(D.12) 
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cr pmn m°'pn = 2 a mn a mn - 9 a , (D.13) 

<? V {ab)q<? Vq = Vpq(ab)V Pq ~ 6 CT(T ab - h>£} <Jpq<J Pq ~ (Jap(T P b (D.14) 

° m e mah h = -30aa a + 9a pq a pqa , (D.15) 

o pq <? pqha h = -24aa a + 7a pq a pqa . (D.16) 



Finally, recall that for any tensor t ab 



V a t b c = D a t b c + - C«Se t b e - - C(% t e c + 1 C( 2 )^ e t b e - 1 C<( 2 ^ t e c + . . . , (D.17) 

where expressions for C' 1 '^ and C^ 2 ' c ab are 

C (1) ^ = -h$ l > c <Ta-hV>°* h + o*h<®, (D.18) 

C^S* = -2 a (4°) c a fe + fcW c a a - /»£) a c ) + I (i? fl ^) c + D^i 2 ) c - D^J) . (D.19) 

These expressions were also given in equations (D.12) and (D.13) of [23]. 
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